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Abstract 

We apply non-extensive methods to the statistical analysis of fully 
developed turbulent flows. Probability density functions of velocity 
differences at distance r obtained by extremizing the Tsallis entropies 
coincide well with what is measured in turbulence experiments. We 
derive a set of relations between the hyperflatness factors F m and the 
non-extensitivity parameter q, which can be used to directly extract 
the function q(r) from experimentally measured structure functions. 
We comment on various non-extensive methods to calculate the mo- 
ment scaling exponents £ m . 



1 Introduction 



Consider a suitable observable in a fully developed hydrodynamic turbulent 
flow. For example, this may be a longitudinal or transverse velocity differ- 
ence, a temperature or a pressure difference. The dynamics is effectively 
described by some highly nonlinear set of model equations — for example the 
Navier-Stokes equation. Since nobody is able to solve this equation exactly, 
it is desirable to find some effective statistical description using methods from 
generalized statistical mechanics. 

The underlying idea is quite similar to what was going on more than a 
hundred years ago, proceeding from classical mechanics to ordinary thermo- 
dynamics. Although nobody was able to 'solve' the classical iV-body problem 
with N > 3 exactly, one still was able to develop quite a successful effective 
theory of a gas of 10 23 particles by extremizing the Boltzmann-Gibbs entropy. 
For fully developed turbulent spatio-temporal chaotic systems, the ordinary 
Boltzmann-Gibbs statistics is not sufficient to describe the (non-Gaussian) 
stationary state. But still we can try to develop an effective probabilistic 
theory using more general information measures. 

The idea to start from an extremization principle for turbulent flows 
is actually not new. For example, Cocke M has presented some work on 
turbulence where he extremizes the Fisher information. Castaing et al. 
also start from another extremum principle to derive probability densities 
in fully developed turbulence. Here we will work within a new approach to 
turbulence ||, f§, [f| based on extremizing the Tsallis entropies f| |7|, || 



The pi are the probabilities of the various microstates of the physical system, 
and q is the non-extensitivity parameter. The ordinary Boltzmann-Gibbs 
entropy is obtained in the limit q — > 1. 

Generally, the Tsallis entropies are known to have several nice properties. 
They are positive, concave, take on their extremum for the uniform distri- 
bution and preserve the Legendre transform structure of thermodynamics. 
On the other hand, they are non-extensive (non-additive for independent 
subsystems) . 

Extremizing S q under suitable norm and energy constraints, one arrives 




(1) 
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at a generalized version of the canonical distribution given by 

Pi = ^(1 + (q-l)(3e i )^, (2) 

where 

Z 9 = £(l-(l-g)/^)^ (3) 

i 

is the partition function, {3 = l/(kT) is a suitable inverse temperature vari- 
able, and the 6j are the energies of the microstates i. Ordinary thermody- 
namics is recovered for q — > 1. One can also work with the escort distri- 
butions ||, defined by Pj = pf/J^Pi- If P is allowed to depend on q, the 
escort distribution is of the same form as eq. (0), with a new q' defined by 
q/(q-l)=:l/(q'-l). 

All we have to decide now is what we should take for the effective en- 
ergy levels 6i in the turbulence application. This depends on the problem 
considered. For example, turbulence in different dimensions ought to yield 
different effective energy levels. Moreover, different observables will also lead 
to different effective energies. In a turbulent 3-dimensional flow for example, 
temperature differences should be described by different effective energies 
than velocity differences. This is clear from the fact that the experimentally 
observed stationary probability distributions (slightly) differ. 

At this stage one has to turn to some sort of model. In the following we 
will consider a simple local model for longitudinal velocity differences that 
seems to reproduce the statistics of true turbulence experiments quite well 
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2 Perturbative approach to chaotically driven 

systems 

The model is based on a generalization of the Langevin equation to deter- 



ministic chaotic driving forces [10|-||12||. We denote the local longitudinal 
velocity difference of two points in the liquid separated by a distance r by u. 
Clearly u relaxes with a certain damping constant 7 and at the same time is 
driven by deterministic chaotic force differences F c h aot {t) in the liquid, which 
are very complicated. Hence a very simple local model is 

u = -ju + F chaot (t). (4) 
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The force F c h ao t{t) is not Gaussian white noise but a complicated determin- 
istic chaotic forcing. It changes on a typical time scale t, which is smaller 
than the relaxation time 7 _1 , so 7r is a small parameter. We may effectively 
discretize in time and consider the rescaled kick force 

oo 

F chao t(t) = (7r) 1/2 E x n6(t - nr), (5) 

n=0 

where the the iterates of an appropriate stroboscopic chaotic map T. 

Integrating eq. (f|) one obtains 

%n+l = T[X n ) 

U n+1 = \u n + y/jTX n+ i, (6) 



where u n := u(nr+0) and A := e 7r . For 7T — > 0, t = nr finite, and so-called 



yj-mixing deterministic maps T it has been shown [ ID | that the w-dynamics 
converges to the Ornstein-Uhlenbeck process, regarding the initial values xq 
as random variables. Hence the invariant density of u becomes Gaussian in 
this limit. For finite 7T, on the other hand, the invariant density is non- 
Gaussian. It can have fractal and singular properties if 7T is large Hl2| . But 



for small 7r it approaches the Gaussian distribution provided T is y9-mixing. 

The route to the Gaussian limit behaviour has been investigated in detail 
in [|T^[ for many different chaotic maps T. Within a well defined universality 
class (defined by square root scaling of the first order correction) it was found 
that for small enough 7T the invariant density of u is always given by 

p(u) = 1 e -KW7^-^ 3 ) + 0( 7 r), (7) 

provided the variable u is rescaled such that the variance of the distribution 
is 1. c is a non-universal constant. This means, not only the Gaussian limit 
distribution is universal (i.e. independent of details of T), but also the way 
the Gaussian is approached if the time scale ratio jt goes to 0. Much more 
details on this can be found in [|TjJ . 

It is now reasonable to assume that also the local chaotic forces acting 
on longitudinal velocity differences in a turbulent flow lie in this universality 
class. We can then use eq. (0) to construct effective energy levels ^ for the 
non-extensive theory. 
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3 Constructing effective energy levels 



Eq. (|7]) corresponds to a Boltzmann factor 

Pi = (8) 

with (3 = l/(kT) = 1, Z = V2tt, and energy q formally given by 

£i=\u 2 ~ c^{u - l -v?) + O ( 7 r) . (9) 

The main contribution is the kinetic energy |w 2 , but in addition there is also 
a small asymmetric term with a universal w-dependence. 

The complicated hydrodynamic interactions and the cascade of energy 
dissipating from larger to smaller levels is now expected to be effectively 



described by a non-extensive theory with the above energy levels (see [|14|] for 
a related cascade model). We obtain from eq. (0) and (|5|) the formula 



1 / „, _ /l o 1 



p(u) = — (l + f3(q-l)(-u 2 - c^n(u - -u") + 0( 7 r) . (10) 



Z q V v ' V2 v ' v 3 
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This equation is in very good agreement with experimentally measured prob- 
ability densities. For detailed comparisons with various turbulence experi- 
ments, see H [15|, |l^, [l7|l . The parameter (3 is determined by the condition 



that the distribution should have variance 1. For 7 r = this is achieved for 
/? = 2/(5-3g). 

Note that we have identified a small parameter 7 r in our approach. It 
is the ratio of two time scales — that of the local forcing and that of the 
relaxation to the stationary state. One may conjecture that it is related to 
the inverse Reynolds number R^ 1 ||. The turbulent statistics is determined 
by a kind of effective non-extensive field theory with the formal coupling 
constant ^J^ff. A perturbative approach is possible since ^fyf is small. In 
fact, for the dynamics ([]) one can work with analogues of Feynman graphs 
related to higher-order correlations of the chaotic dynamics |18| , |19|| . Eq. (|10|) 
is just obtained by first-order perturbation theory — the complete theory is 
the infinite-order theory taking into account all orders of \Ht- But first we 
have to understand the 'free' turbulent field theory obtained for 7 r = 0. Here 
almost everything can be calculated analytically. 



4 



4 The 'free' turbulent field theory 

If 7T = the moments can be easily evaluated. In || we obtained 

(H1 = /, f MH"*=( T ) ^fprrf (») 

with fc defined as := l/(q — 1) (A; needs not to be integer). The beta 
function is defined as 

™-WSr (12) 

We now show that this formula for the moments can be significantly simpli- 
fied. First, we obtain from the definition of the beta function 

<H"> = ( 7 ) 1 <13) 

Generally, one has 

n— 1 

r(x + n) = ro) J] (s + j) (14) 

i=o 

for natural numbers n. Suppose y =: n e N, then it follows from eq. (|i~3|) 
and (0) 

V/3/ ,i 1 2j + 2*-m-l 



In particular, we obtain for the first few moments 



<" 2 > = (f)sb (16) 

^ (f) 2 2^5 2^3 W 

(19) 
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This can be used to evaluate the complete set of hyperflatness factors F, 
defined as 

(u 2m ) 

W) 



We obtain 



and generally 



F m = Swf (20) 



Fx = 1 (21) 
2/fc — 3 

2k^5 (22) 

* = 3 ' 5 (^W^) 3 (23) 

F ' = 3 ' 5 ' 7 (2fc-9) ( (2fc-7)(2A;-5) (24) 



(2A;-3) m -' 
1?= 5 +1 (2^ 

(j odd). Note that all hyperflatness factors are independent of f3. 



F -° (2m -' )!! igPM (25) 



5 Extracting g(r) from experimentally mea- 
sured structure functions 

The great advantage of the hyperflatness factors is that they yield a direct 
way to estimate the r-dependent non-extensitivity parameter q(r) from ex- 
perimentally measured structure functions (u m )(r). Eq. (p2|) yields 



or 



equivalent to 



(see also 



F _6fc-9_ 15-9g 



7^2 - 15 



Another relation follows from eq. (f23|). 



2F 3 - 30 

equivalent to 



6F 3 - 45 ± JF$ + I2OF3 , (29) 



8F 3 - 75 ± JF? + I2OF3 

q = — V , \ (30) 

6F3 - 45 ± JF? + I2OF3 



In fact, each hyperflatness factor F m with m > 2 yields a relation for (or 
q), and all relations are the same in case the 'free' turbulence theory is exact. 

Given some experimentally measured hyperflatness structure functions 
F m { r ) one can now determine the corresponding curves q(r). The less these 
curves differ for the various m, the more precise is the zeroth-order (free) 
turbulence theory. For examples of such experimentally measured curves 
q(r), see Generally, the hyperflatness factors are complicated functions 



of both the Reynolds number R\ and the separation distance r, and so is 
q(R\,r). But one may conjecture that for R\ — > 00 one obtains a universal 
function q*(r). How it looks like in the entire r-range is still an open question. 



6 Some remarks on the scaling exponents ( m 

The precise values of the scaling exponents ( m , which describe the scaling 
behaviour of the structure functions (u m ) ~ r^ m in the inertial range, are still 
a rather controversial topic in turbulence theory. In fact, it is not at all clear 
whether there is exact scaling at all or just approximate scaling, whether the 
higher moments exist at all or not, how reliable the experimentally measured 
higher-order exponents are, and what the effects of finite Reynolds numbers 
are. On the theoretical side, a variety of models and theories have been 
suggested (see, e.g., [^, [Hj]) but a true breakthrough convincing a significant 
majority of scientists working in the field seems not to have been reached at 
the moment. 

To 'derive' values for the scaling exponents using methods from non- 
extensive statistical mechanics, one needs additional assumptions — just as 
in all other models and theories dealing with the scaling exponents. In § a 
logarithmic depence of k on r was suggested to derive the following formula 
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for the scaling exponents 



Here fc = 1/(5 — 1) is the average value of the non-extensitivity parameter in 
the inertial range. On the other hand, Kolmogorov's lognormal model (the 
K62-theory f2~2| ) predicts 



m ( n\ 1 2 



where /1 is the intermittency parameter. Now, asuming that k is rather large 
we can expand the logarithms in eq. fl3"Tp. The linear terms cancel, and the 
first non-trivial terms are the quadratic ones. Neglecting the higher-order 
cubic terms, one precisely obtains from eq. (PH ) the result Q3"2]), identifying 

" = ln2( 2 Ll)^ (33) 

It is encouraging that the simplest non-extensive model assumptions lead to 
an old theory that has a long tradition in turbulence theory (though it is 
known that K62 cannot be correct — it is, however, a good approximation 
for m not too large). The value of the intermittency parameter comes out 
with the correct order of magnitude (/i »i 0.2) if one chooses k rs 6, and the 
same k fits the experimentally measured probability densities correctly. 

On the other hand, starting in the derivation leading to eq. ( |3~TD from 
a different dependence of k on r rather than the logarithmic one, one also 
obtains different formulas for the scaling exponents. So the final answer to 
the question of the scaling exponents is still open. 



Actually, in || eq. (|31| ) was derived using the 'free' turbulence theory, but 
clearly the correct theory is the infinite-order theory. As one can easily see, 
the infinite-order theory yields probability densities p(u) living on a compact 
support for arbitrary small but finite jt. Indeed, iterating eq. @ one obtains 
u n = X u uq + ^/ttZ)"=i X n ~^Xj. Hence, for any chaotic dynamics T bounded 
on some finite phase space (say [-1,1]) one obtains for n — > 00 from the 
geometric series the rigorous bound \u\ < ypfr/(l — A) ~ \jJyf. Thus 



for the infinite-order theory all moments exist for arbitrary small but finite 
7T — in contrast to the 0-th order theory, where only moments (u m ) with 
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m < 2k — 1 exist, since here the densities live on a non-compact support and 
there is polynomial decay for large \u\. This once again shows how delicate 
the problem of the existence of the moments and of the scaling exponents 
in general is. For densities living on a compact support, as provided by the 
infinite-order theory, one expects a linear asymptotics of ( m for large m. 

Generally, the way non-extensive statistics can be used to 'derive' the 
scaling exponents ( m is not unique. An alternative approach, based on an 
extension of the multifractal model and q- values smaller than 1, has been 
suggested by T. and N. Arimitsu [f|. In their model asymptotically there is 
a logarithmic dependence og ( m on m. 
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